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Absolute total cross sections for electron capture between slow, highly charged ions and alkali 
targets have been recently measured. It is found that these cross sections follow a scaling law 
with the projectile charge which is different from the one previously proposed basing on a classical 
over-barrier model (OBM) and verified using rare gases and molecules as targets. In this paper 
we develop a "semi-classical" (i.e. including some quantal features) OBM attempting to recover 
experimental results. The method is then applied to ion-hydrogen collisions and compared with the 
result of a sophisticated quantum-mechanical calculation. In the former case the accordance is very 
good, while in the latter one no so satisfactory results are found. A qualitative explanation for the 
discrepancies is attempted. 

PACS numbers: 34.70+e, 34.10. +x 



I. INTRODUCTION 

The electron capture processes in collisions of slow, highly charged ions with neutral atoms and 
molecules are of great importance not only in basic atomic physics but also in applied fields such as 
fusion plasmas and astrophysics. 

In the past years a number of measurements have been carried on the collisions between highly 
charged ions and rare gases Q or molecules ||, in which one or several electrons were transferred 
from the neutral target to a charged projectile: 

A +q + B -> A {q - J)+ + B J+ . (1) 

Their results-together with those from a number of other laboratories-yielded a curve which can 
be fitted within a single scaling law (a linear relationship) when plotting cross section a versus 
projectile charge q: it is almost independent of the projectile species and of the impact velocity 
v (at least in the low-speed range v < 1 au). When one extends experiments to different target 
species, the same linear relation holds between a and q/It, with It the ionization potential of the 
target |§. 

It is found that this scaling law could to be predicted, in the limit of very high projectile charge, by 
a modification of an extended classical over-barrier model (ECBM), allowing for multiple electron 
capture, proposed by Niehaus Quite recently a confirmation of this scaling has come from a 
sophisticated quantum-mechanical calculation 

Similar experiments were carried on more recently for collisions between ions and alkali atoms 
Q. The results show that the linear trend is roughly satisfied, but the slope of the straight line is 
grossly overestimated by the ECBM: in Fig. 1 we show some data points (stars with error bars) 
together with the analytical curve from the ECBM (dashed curve) which, for one-electron atoms, 
is written J1B| 

a = 2.6 x W 3 q/lf [l(T 20 m 2 ] (2) 

(It in eV). It should be noticed that experimental data are instead well fitted by the results of a 
Classical Trajectory MonteCarlo (CTMC) code Q. 

The ECBM of ref. j^] works in a simplified one-dimensional geometry where the only physically 
meaningful spatial dimension is along the internuclear axis. It does not take into account the fact 
that the electrons move in a three-dimensional space. This means that only a fraction of the electrons 
actually can fulfil the conditions dictated by the model. For rare gases and molecules, which have a 
large number of active electrons, this can be not a trouble (i.e., there are nearly always one or more 
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electrons which can participate to the collision). For alkali atoms with only one active electron, on 
the other hand, an overestimate of the capture probability by OBM's could be foreseen. 

With present-days supercomputers there are relatively few difficulties in computing cross sec- 
tions from numerical integration of the time-dependent Schrodinger equation (e.g. refer to ref. Q). 
Notwithstanding this, simple models are still valuable since they allow to get analytical estimates 
which are easy to adapt to particular cases, and give physical insight on the features of the problem. 
For this reason new models are being still developed HM. 

In this paper we present a modified OBM which allows to get a better agreement with the experi- 
mental data of ref. |7|. 



II. THE MODEL 

We start from the same approach as Ostrovsky |j| (see also jj^]): be r the electron vector relative 
to the neutral atom (T) and R the internuclear vector between T and the projectile P (see Fig. 
2 for a picture of the geometry: it is an adaptation of Figure 1 from ref. |]]). Let us consider the 
plane containing the electron, P and T, and use cylindrical polar coordinates (p, z, (j>) to describe the 
position of the electron within this plane. We can choose the angle cj> — and the z direction along 
the internuclear axis. We will assume that the target atom can be described as an hydrogenlike 
atom, which is not a bad approximation when dealing with alkali atoms. 
The total energy of the electron is 

E=£ + U= P 1 Z 7 
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vV + z 2 ^J P 2 + {R- z y 

Z p and Zt are the charge of the projectile and the effective charge of the target seen by the electron, 
respectively, and we are using atomic units. 
We can also approximate E as 



E(R) = -En 



R 



E n is given by the quantum-mechanical value: E n = Z?/(2n 2 ). This expression is asimptotically 
correct as R — > oo. 

On the plane (e, P, T) we can draw a section of the equipotential surface 

U(z,p,R) = -E n -^- . 

This represents the limit of the region classically allowed to the electron. When R — > oo this region is 
decomposed into two disconnected circles centered around each of the two nuclei. Initial conditions 
determine which of the two regions actually the electron lives in. 

As R diminishes there can be eventually a time where the two regions become connected. It is easy 
to solve eq. (|B|) for R by imposing that p m = and that there must be an unique solution for z 
with < z < R: 



Zt + 2yJ ZtZp 

Rm = p ■ 

-Em 

In the spirit of OBMs it is the opening of the equipotential curve between P and T which leads 
to a leakage of electrons from one nucleus to another, and therefore to charge exchange. Along the 
internuclear axis the potential U has a maximum at 

R ^ 

z — zo = R- 



Whether the electron crosses this potential barrier depends upon its initial conditions. These are 
chosen from a statistical ensemble, which we will leave unspecified for the moment. Let Nn be the 
fraction of trajectories which lead to electron loss at the time t and W(t) the probability for the 
electron to be still bound to the target, always at time t. The fraction of losses in the interval t,t + dt 
is given by 

dW{t) = -Nn^W{t) , 

J- em 



with T em the period of the electron motion along its orbit. A simple integration yields the leakage 
probability 



P ; = l-exp^-^- j Nndtj . (9) 

In order to actually integrate Eq. (^) we need to know the collision trajectory; an unperturbed 
straight line with b impact parameter is assumed: 



R = \A 2 + (vt) 2 . (io) 

At this point it is necessary to give an explicit expression for Nn ■ The electron is supposed to be in 
the ground state (n = 1,1 = m = 0). T em becomes therefore |n[ 

T em = 2tt/Z? . (11) 

Ref. |]] adopts a geometrical reasoning: the classical electron trajectories, with zero angular momen- 
tum, are ellipses squeezed onto the target nucleus. The only trajectories which are allowed to escape 
are those whose aphelia are directed towards the opening within the angle ±# m . The integration 
over this angle yields an analytical expression for Nn (Eq. 17 of ref. In Fig. 1 we show the 

results obtained using Ostrovsky's model ( dotted curve-eqns. 8,17 of ref. 0) Notice that from 
direct inspection of the analytical formula, one sees that the scaling law is not exactly satisfied, at 
least at small values of the parameter q/It, and this is clearly visible in the plot. The result is 
almost equal to the scaling (□). 

The present approach is based on the electron position instead than on electron direction . The 
recipe used here is (I) to neglect the dependence from the angle: all electrons have the same proba- 
bility of escaping, regardless of their initial phase. Instead, (II) the lost electrons are precisely those 
which, when picked up from the statistical ensemble, are found farther from nucleus T than the 
distance zo- 

poo 

N n = f(r)dr , (12) 

J ZQ 

with f(r) the electron distribution function. 

There is not a unique choice for f(r): the (phase-space) microcanonical distribution 

/( r ,p)o(ik + ^-^ (13) 

(5 is the Dirac delta) has been often used in literature since the works [^j as it is known that, 
when integrated over spatial coordinates, it reproduces the correct quantum-mechanical momentum 
distribution function for the case of the electron in the ground state [[l3j (more recently the same 
formalism has been extended to Rydberg atom s Hi ). After integration over momentum variables 
one gets instead a spatial distribution function J15|| 



/mc(r) = ^ 2Z *) 3/ W l-gi, r <2/Z t (14) 
7T V r 2 

and zero elsewhere (The lowerscript "mc" is to emphasize that it is obtained from the microcanonical 
distribution). However, this choice was found to give poor results. It could be expected on the basis 
of the fact that (^) does not extend beyond r = 2/Z t and misses therefore all the large impact- 
parameter collisions. In the spirit of the present approach, it should be instead important to have an 
accurate representation of the spatial distribution. We use therefore for f(r) the quantum mechanical 
formula for an electron in the ground state: 

fu(r) = 4Z?r 2 exp {-2Z t r) (15) 

which, when substituted in (fl^), gives 



1 Beware of a small difference in notation between the present paper and j8|: here we use an effective charge for the target, 
Zt = \f2E n , while |8j uses an effective quantum number nt = l/y/2E n with the effective charge of the target set to 1. 
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N n = [1 + 2z Z t + 2{z Z t ) 2 ] exp (-2z Z t ) . 

Since the choice for f(r) does not derive from any classical consideration, we call this method a 
"semi-classical" OBM. 

Notice that, in principle, one could go further and compute f(r) from a very accurate wavefunction, 
fruit of quantum mechanical computations (see (lfi)), but this is beyond the purpose of the present 
paper (it could be worthy mentioning a number of other attempts of building stationary distributions 
f(r), mainly in connections with CTMC studies, see [fl7l-|is|l). 

The f(r) of Eq. (|l^) does not reproduce the correct momentum distribution, nor the correct energy 
distribution (which could be obtained only by using eq. jl3|). However, it is shown in jl5| that this 
choice gives an energy distribution for the electrons, f(E), peaked around the correct value E n , and 
< E >= E n , where < . . . > is the average over f{E). 

Some important remarks are to be done here. First of all, a question to be answered is: why use an 
unperturbed distribution, when the correct one should be sensitively modified by the approaching of 
the projectile. The answer is, obviously, that this choice allows to perform calculations analitically. 
We are doing here a sort of classical counterpart of a quantum-mechanical Born calculation: there, 
too, the matrix elements are computed as scalar products over unperturbed states, regardless of 
any perturbation induced by the projectile. In the following, however, some considerations about 
possible improvements over this simple approximation will be done. 

A second question regards the meaning of the factor dt/T em in eq. (^): in Ostrovsky's paper 
this is the fraction of electrons which enter the loss zone during the time interval dt and is valid 
under the hypothesis of a uniform distribution of initial phases of the electrons. In our case this 
this assumption ceases to be valid: electrons actually spend different fractions of their time at 
different radial distances from T, depending on their energy. We will do a (hopefully not too severe) 
assumption by assuming that, on the average, the expression (0) still holds. 



III. RESULTS 
A. Iodine - Cesium 

This study has been prompted by the ion-atom experiments of j?| : first of all, therefore, we apply 
the above model to the process of electron capture 

l q+ + Cs^ + Cs+ 

with q = 6 -7- 30. Impact energy is 1.5 x q keV [Q. The ionization potential of Cesium is I t = 3.9 
eV. Solid line in Fig. 1 is the result of the present model: the agreement is fairly good. 



B. Bare ions - Hydrogen 

As second test, we have computed cross section for captures 

H + O s+ — > H + + 7+ 

and 

H + He 2+ -> H+ + He+ 

and compared it with similar calculations done using the molecular approach by Harel et al pofl . 
The results are summarized in fig. 3. There is a sharp discrepancy in the behaviour for v — > 0, 
where the present model predicts an increasing cross section. At very low speed it is the concept 
itself of atomic distribution function which becomes questionable, and molecular aspects become 
important. Besides, quantum effects such as the discreteness of the energy levels also play a major 
role and are completely missed by this approach. In the higher velocity part, the present model 
underestimates the more accurate value by a factor 2 for process (^), but the error is much less, 
just 25 %, for process ([if]). These two ions have been chosen ad hoc: they correspond to values of 
the ratio Zt/Z p = 1/8 and 1/2 respectively. In the (I, Cs) test this ratio ranged from « 1/12 to 
~ 1/60 depending upon the projectile charge. This means that in the former case the perturbation 
of the projectile on the electron distribution function is comparable to the (I, Cs) case, while in the 
latter it is much less. We expect the electron distribution function to be more and more perturbed 
as Zt/Z p -* 0. 



IV. SUMMARY AND CONCLUSIONS 



We have developed in this paper a very simple OBM for charge exchange. It exploits some features 
of the quantum mechanical version of the problem, thus differing from similar models which are 
solely classical. The agreement with experiment is much better than previous calculations where a 
comparison could be made. It is far from excellent, but reasons for the (partial) failure have been 
suggested. 

As it stands, the model is well suited for one-optical-electron atoms (since it uses hydrogen-like 
wavefunctions) , therefore we do expect that other classical OBM's can still work better in the 
many-electrons targets studied in previous experiments. 

Some improvements are likely to be added to the present model: a possible line of investigation 
could be coupling the present method with a very simplified calculation of the evolution of the 
wavefunction, using quantum mechanics. From this one should not compute the / as coming from 
a single state, but as a linear combination including also excited wavefunctions (the relative weights 
in the combination should be given by the quantum mechanical calculation) . Work in this direction 
is currently underway. 
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FIGURE CAPTIONS 
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FIG. 1. Comparison between experimental data and prediction from models for electron capture cross section of process 
([[?]). Stars, experiment with 20% error bar; dashed line, scaling law from Niehaus (Eq. ^); dotted line, Ostrovsky's scaling 
law; solid line, scaling law from present model, a is in l(T 20 m 2 , I t in eV. 
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FIG. 2. Geometry of the scattering. P and T are the projectile and target nucleus respectively. The enveloping curve shows 
a section of the equipotential surface U = E, i.e., it is the border of the region classically accessible to the electron. R is the 
internuclear distance. The parameter p m is the radius of the opening which joins the potential wells, 6 m the opening angle 
from T; zq is the position of the potential's saddle point. 
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FIG. 3. Capture cross section versus impact velocity. Upper, H-0 8+ collisions; lower, H-He 2+ collisions. Diamonds, data 
from ref. 20; solid line, present model. 
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